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PREFACE 




This document serves as a summary of the NASA Lewis Research Center (LeRC) Workshop 
on Forced Response in Turbomachinery in August of 1993. The workshop was sponsored by the 
following NASA organizations: Structures, Space Propulsion Technology, and Propulsion 
Systems Divisions of the NASA Lewis Research Center and the Aeronautics and Advanced 
Concepts & Technology Offices of NASA Headquarters. In addition, the workshop was held in 
conjunction with the GUIde (Govemment/Industry/Universities) Consortium on Forced 
Response. The workshop was specifically designed to receive suggestions and comments from 
industry on current research at NASA LeRC in the area of forced vibratory response of 
turbomachinery blades which includes both computational and experimental approaches. There 
were eight presentations and a code demonstration. Major areas of research included aeroelastic 
response, steady and unsteady fluid dynamics, mistuning, and corresponding experimental work. 

Forced response in turbomachinery blades is one of the most important issues in blade 
design. Controlling the forced response of turbomachinery blades is crucial for reliable operation 
of rocket and aircraft engines. Yet a major vacuum exists in predicting dynamic stresses in 
turbomachinery blades due to forced response. Furthermore, with the next generation engines 
being conceived for advanced aircraft, such as the high speed civil transport, forced response 
problems will have a major impact on their design and development. The acuteness of the problem 
is further emphasized by the willingness of all the engine companies to work together in the 
GUIde Consortium. 

Representatives from the U.S. aerospace engine companies (Textron Lycoming, General 
Electric, Rocketdyne, Pratt & Whitney Rocket Engine Group, Pratt & Whitney Aircraft Engine 
Group, and United Technologies Research Center), numerous universities, the U.S. Air Force, and 
NASA LeRC attended the one-day workshop. In addition, the representatives from these 
companies, government agencies, and universities form the nucleus of the GUIde Consortium. 
The consortium was created to perform basic research on forced response of turbomachinery. This 
research work will improve the durability of turbomachinery and prevent engine failures. 

An excellent dialogue took place in which many suggestions were given in a round-table 
discussion following each presentation. Action items from the participants were evaluated and 
incorporated into future work. An industry sub-committee was created to guide the FREPS 
(Forced REsponse Prediction System) development. Currently, FREPS is being modified to 
formulate a workstation version in order for industry to access the code for their design 
requirements. The industry representatives unanimously decided to participate in the Forced 
Response workshop every two years to evaluate NASA’s research efforts in this area. 
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PROGRAM OVERVIEW 


G.L. Stefko 

NASA Lewis Research Center 
Cleveland, Ohio 44135 

PROGRAM OVERVIEW 

o FREPS SYSTEM 
o EARTH-TO-ORBIT PROGRAM 
o NASA LEWIS TEAM 
o EXTERNAL TEAM MEMBERS 
o FREPS PLANNED DEVELOPMENT 
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Blade Forced REsponse Prediction System 
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EARTH-TO-ORBIT PROPULSION 

- TURBOMACHINERY SUBSYSTEMS - 
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DESIGN TOOLS 
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FREPS—OBJECTIVES AND OVERVIEW 


D.V. Murthy* 

NASA Lewis Research Center 
Cleveland, Ohio 44135 


Forced Vibratory Response in Turbomachinery 


Consequences 


> Decreases the Fatigue Life due to 
HCF Failures 

0 Increases Development and 
Maintenance Cost 

£> Imposes Operational Restrictions 


Future Trends Point to 
Growing Problems 


P3P Higher Power-to-Weight Ratio Goals 

Lower Damping Designs (e.g. Blisks) 

H5F Increased Exposure to 

Resonances (e.g. Low AR Blades) 


Rocket Engines vs Air Breathing Engines 
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Turbine Inlet Temp, deg f 

1540 

2600 

Blades Cooled? 

No 

Yes 

Blade Metal Temp, deg F 

1500 

1500 Cruise 

Rotor Tip Speed, fps 

1850 

1650 Take off 

HP/Rotor Blade 

500 

600 

Material 

Super Alloy - MAR-M 24 

6 Super Alloy - MAR-M 246 

Useful Life, hours 

0.1 to 7.5 

8000 


NASA Resident Research Associate at Lewis Research Center. 
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Origins of Unsteady Aerodynamic Excitation! 



Inlet Flow 
Distortion and 
Turbulence 


* 


Compressor 
Surge and 
Rotating Stall 


Origins of Unsteady Mechanical Excitation | 



Blade Tip-Casing 
Contact 



Disk Flexibility 


* 

* 



* 



^ Shaft and Gear 
Meshing 


^ Foreign Object 
Damage 


Cooling Jets 
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FREQUENCY (KHZ) 


| Campbell Diagram] 
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FREPS Overcomes the Above Limitations and Provides an 
Efficient Tool for Design Environment 



Forced REsponse Prediction System 


Objective 

To Develop a Design Tool for Predicting the Forced Vibratory 
Response of Turbomachinery Blades to Unsteady Excitations 

To Validate using Test/Experimental Data 
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HIM I 


Equations of Motion 



Mechanical/Material Damping External Forcing Function from 

Forces from Damping Model Mechanical Excitation Model 

/W, 



Elastic and Inertia Forces Unsteady Aerodynamic Forces 
from Structural Model from Aerodynamic Model 




ROTOR 

BLADE 

AERODYNAMIC 

DAMPING 

Tuned 

Typical Section 

Fully Linear 

Structural 

identical Blades 

Lumped Parameter 

Empirical 



Linearized 

Mistuned 

Beam 

Potential 

Aerodynamic 

Blade-to-Blade 

Id variation 

Euler 

Variations 

Shell/Plate 

Friction 

Nonlinear 

Rigid Disk 

2d variation 

Potential 

Euler 

impact 

Flexible Disk 

Finite Element 
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Goals | 

Accurate Predictions 
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Detailed Models 

Steady Flow 
Unsteady Flow 


Rotor 

Blade 


| Rotor Structural Model: Tuned or Mistimed ?] 

Tuned Model is More Suitable to Design Procedures: 

Model One Blade Only 
<y" Smaller Problem Size 

Usually Designer's Preference 

Mistuned Model Difficult to Incorporate in Design Procedures: 

Si Much Larger Problem Size 
Si Greater Complexity 

Si Mistuning Pattern / Strength Unknown at Design Time 
Si Mistuning Pattern / Strength Vary from Rotor to Rotor 
Si Monte Carlo Simulations are Very Expensive 


9 





Blade Structural Model 



Aerodynamic Model 
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Computational Efficiency 


Solution Methods 

Linearized Dynamic Analysis 
Modal Solution 

Statistical Treatment of Mistuning 
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Computational Efficiency 


Solution Methods 


Linearized Dynamic Analysis 


Modal Solution 


Statistical Treatment of Mistuning 
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Computational Efficiency 

Solution Methods 

Linearized Dynamic Analysis 

Modal Solution 

Statistical Treatment of Mistiming 





Computational Efficiency 


Solution Methods 

Linearized Dynamic Analysis 


Modal Solution 



Statistical Treatment of Mistuning 
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UNSTEADY AERODYNAMIC ANALYSES FOR TURBOMACHINERY 


AEROELASTIC PREDICTIONS 

J.M. Verdon, M. Barnett, and T.C. Ayer 
United Technologies Research Center 
East Hartford, Connecticut 06108 
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UNSTEADY AERODYNAMIC 
ANALYSES 


• Applications 

- Aeroelastic: blade flutter and forced vibration 

- Aeroacoustic: noise generation 
— Vibration and noise control 

- Effects of unsteadiness on performance 

• Requirements 

- Accuracy /efficiency 

* Realistic operating conditions 

* Arbitrary modes of unsteady excitation 

• Approaches 

— Numerical simulation/analytical modeling 


ASSUMPTIONS 


• Turbulence and transition can be modeled 

=► Reynolds averaged, Navier-Stokes equations 

• High Reynolds number, “attached” flow 

=► Thin-layer Navier-Stokes equations, or 
Inviscid/viscid interaction analyses 


• Small-amplitude unsteady excitations 

=*• Nonlinear steady + linearized unsteady analyses 

• Re —► oo ==> inviscid flow 

- Potential steady background flow => LINFLO 

- Uniform steady background flow => CLT 
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CONTRACT NAS3-25425 


NASA Program Managers: J. Gauntner, G. Stefko 

• Linearized inviscid unsteady aerodynamic analysis: LINFLO 

• Unsteady viscous layer analysis: UNSVIS 

• Steady, inviscid/viscid interaction analysis: SFLOW-IVI 

• Coupled SFLOW-IVT/LINFLO analysis 


EXAMPLE 

CONFIGURATION 
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UNSTEADY EXCITATIONS 


S_oo(x, t) 







40 


PI- oo(x, t) *j 

\m 

PL+oo(x,t) 

Y 

! 

<-co(x,0 

W-t 





• Far-field conditions (uniform mean flow) 

*(x, 0 = Pe{s_ ro expf^K-oo • x + wt)]}, £ < £_ 

<(x, t ) = Re{<_ 0O expl^/c-oo • x + wt)]}, { < £_ 

t) = Re{p/,Too exp[—A foo£ + '(ktoo * x + wt)]}, f > & 


LINEARIZED INVISCID ANALYSES 


• Linearization 

P(x, t) = P(x) + Re{p(x) exp(turt)} + ... 


• Nonlinear BVP for steady background flow 

• Linear variable-coefficient problem for each Fourier 
component of first-order unsteady flow 

- Time independent 

- Surface conditions imposed at mean surfaces 

- Analytic far-field solutions for s, C, and P 

- Single extended blade-passage solution domain 

P(x + mG,t) = P(x) + Re{p(x) exp[»(u;< + mu)]} + ... 

• Prescribed quantities: 

r B> ^—oo C—oo> an d Pirfoo 
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LINFLO 


• Unsteady perturbation of a potential mean flow 

• Steady flow: V • (pV$) = 0 

• Unsteady velocity decomposition: v = V(<^ + ^.) + v R 

— p = -pD<t>/Dt 

— v ■ v/j = 0 far upstream 

I 

— D<t>,/Dt = 0; (V^. + v R ) ■ n = 0 on B m k W m 

• Entropy & rotational velocity: X = Aer+^e N -♦ x as £ -* -oo 

s(x) = s_oo exp(«/c_oo ■ X) 

V«(x) = [V(X • A_oo) + s-«,V$/2] x exp(iK_oo • X) 

• Unsteady velocity potential 

D(A~ 2 D<t>/Dt)/Dt -p~ l V- (pV<f>) = p~ l V- \pV4>.] 
where <j>, = F(A-oo, 'S>) exp(i/e_ 00 • X) 

• Surface conditions: 

- Blades: V<£ • n = /(r b) 

- Wakes: [D4>/Dt] = 0 and [V^] - n = 0 

- Shocks: [pV<f> + pV$>] • n = /(rsh • n, V$); rsh • n = -[^]/|$ n J 

• Far field conditions: 

- 4>i &oo prescribed; 4>r,*o o must be determined 

- Analytic far-field solutions for <j> = <f>i + <j> R 
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NUMERICAL SOLUTION DOMAIN 


• Extended blade-passage region of finite extent in axial- 
flow direction 



NUMERICAL APPROXIMATION 


• Implicit, least-squares, finite-difference model 

{C<t>)o « (L<f>) 0 = q°<j>o + J2 AnOm - 4>o) 

m=l 

• Transonic differencing strategies 

• Cascade, local and composite mesh solutions 

• Direct solution procedure 

- Block tridiagonal system of algebraic equations for 
subsonic flow 

- Block pentadiagonal system for transonic flow with 
fitted shocks 
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AERODYNAMIC RESPONSE AT A 
BLADE SURFACE 


• Surface pressure (transonic flow): 

P( T B, t) = P{r b) + Re{p 5 (r B ) exp(iwl)} + <) + ••■ 

n 

• Blade motion: r B (x) = E^R,(x) 

i= 1 

• Unsteady airloads: 

9i = 0 qi,rdr = -(f) [P^TT- X e, + Pan - X! r Sh„[-P]n] • R,dr 

• Work per cycle/pressure-displacement function 

W C =0 = ^ w ( r )d T = *0 kn {6’q iT }dT = 7rlm{ £ £*9<} 


EXAMPLE RESPONSE PREDICTIONS 


• Compressor exit guide vane (EGV): 0 = 15 deg, G = 0.6 

- Thick, highly-cambered NACA 0012 airfoils 

- Subsonic flow: = 0.3, = 40 deg 

- Vortical excitation: w = 10, cr = -2tt 

— Acoustic excitation from downstream: 

U — 10, (7 = 0 

• High speed compressor cascade: 0 = 45 deg, G = 1 

- Cambered NACA 0006 airfoils 

- Subsonic flow: = 0.7, = 58deg 

- Transonic flow: M_« = 0.8, = 55 deg 

- SDOF blade motions: <5, = (1,0), u> = 1 

• Lin ear/nonlinear result comparisons 

- NGUST analysis (Navier-Stokes) 

- NPHASE analysis (Euler) 
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COMPRESSOR EXIT GUIDE VANE 

NGUST Computational Grid 



COMPRESSOR EXIT GUIDE VANE 


Steady surface pressure coefficient 


- CASPOF (Pull potential) 

-NGUST (Euler) , 

-NGUST (Navier-Stoke^^-^rrr^r^ 















I VORTICITY WAVE IN AN EGV CASCADE 
Unsteady vorticity, vr-oo = (0.05$, 0), a = -2?r, u = 10.0 



[VORTICITY WAVE IN AN EGV CASCADE 
Unsteady pressure, v R = (0.05?, 0), a = -2tt, w = 10.0 
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VORTICITY WAVE IN AN EGV CASCADE 


FIRST HARMONIC UNSTEADY PRESSURE DIFFERENCE 
vr_ oo = (0.05g,0), (7 = -27T, w = 10.0 



VORTICITY WAVE IN AN EGV CASCADE 


FIRST HARMONIC UNSTEADY PRESSURE DIFFERENCE 


Vr-co = (0.05<j, 0), <? = -27r, uj = 10.0 

- UR .-00 = 0.05 a 

- Vr-co = Q.lOq 

- Vr-co = 0.25g 

- 4-oo = 0.50g 

VISCOUS SIMULATIONS 
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COMPRESSOR EXIT GUIDE VANE 


Unsteady Pressure Response 
p+oo = (0.04,0), lo = 10.0, c — 0.0 

Linearized Invisdd (LINFLO) Navier-Stokes (NGUST) 
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EXIT ACOUSTIC WAVE IN AN EGV CASCADE 


FIRST HARMONIC UNSTEADY PRESSURE DIFFERENCE 
P/,+00 = (0.04P, 0), a = 0, w = 10.0 
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EXIT ACOUSTIC WAVE IN AN EGV CASCADE 


FIRST HARMONIC UNSTEADY PRESSURE DIFFERENCE 
P/,+oo = (0.04P, 0), <r = 0, u = 10.0 


' P/,+CC — 

-P/,+co = 

-P/,+00 = 


0.04P 

0.12P 

0.20P 



NACA 0006 CASCADE 

NPHASE Computational Grid 
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HIGH SPEED COMPRESSOR CASCADE 

Surface Mach Number Distributions 
-Potential,-Euler 



HIGH SPEED COMPRESSOR CASCADE 

Pressure Displacement Function Distributions for 
Torsional Blade Vibrations at a = 2 deg, w = 1 
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HIGH SPEED COMPRESSOR CASCADE 

Work per Cycle versus Interblade Phase Angle for 
Torsional Blade Vibrations at a = 2 deg, u — 1 


□ Linearized Analysis (LINFLO) 

-Nonlinear Euler Analysis (NPHASE) 

Subsonic Flow Transonic Flow 



INVISCID/VISCID INTERACTION 
ANALYSES 

• High Reynolds Number Flow 

• Inviscid region: Euler or potential flow equations 

- Surface conditions modified to account for viscous 
displacement effects 

• Viscous region: Prandtl’s equations 

- Direct solution: P —* 6 

- Inverse solution: 6 —► P 

• Inviscid viscid interaction law 

- Weak interaction =>• sequential solution, pressure 
determined by inviscid flow 

- Strong interaction => simultaneous solution, pressure 
determined by inviscid and viscous flows 
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CASCADE FLOW WITH LOCAL 
REGIONS OF STRONG INTERACTION 



Leading edge 
separation 


M<1 /M>1 / M<1 



Strong- 

interaction 

regions 


V-co 


Shock/boundary-layer 
interaction 


Trailing-edge 

separation 


Voo 


Weak interactions 


SFLOW-IVI: INVISCID REGION 


• Field equation 

- pV<t> = 0 or A 2 V 2 <t> = V$ ■ V(V$) 2 /2 

• Surface b.c.’s account for viscous displacement effects; 
i.e., 

- Blades: V$ • n| 5 = p~ l d(p e u e 6)/ds 

- Wakes: [V<t>] ■ n + | w = p~ l d{p t u t 6 w )/ds 

• Inlet flow conditions prescribed 

• Exit flow conditions determined by Kutta cond. &z 
global mass conservation 
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SFLOW-IVI: VISCOUS REGION 


• Classical Viscous-Layer Eqs. (Boundary layers & Wakes) 

— Weak interaction: 

specify dujds —♦ calc. 6* (direct) 

— Strong interaction: 

specify m = p t u t 8* —♦ calc. u t (inverse) 

• Turbulence and transition 

— Algebraic eddy-viscosity model 

* Blade: Cebeci-Smith w/separation modification 

* Wake: Chang, et al 

— Instantaneous transition 

• Solutions in terms of Levy Lees variables 


SFLOW-IVI: COUPLING 
PROCEDURE 
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NUMERICAL RESULTS 


• Two Cascade Configurations 

- Compressor exit guide vane (EGV) 

— High-speed compressor cascade 

• Effect of Varying Re 

• Comparison with Navier-Stokes solutions 

• Incidence Angle Study (EGV) 


COMPRESSOR EXIT GUIDE VANE 


Effect of Varying Re: - Re = 10 5 

- Re = 10 6 

-Inviscid 



0.100 

0.075 

6 

0.050 

0.025 


0.0 

0.0 0.5 1.0 1.5 2.0 


surface 

wake 

/ 

1 

suction J 
surface / 

\ /; 
pressure /V 
surface /: 
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COMPRESSOR EXIT GUIDE VANE 


Comparison with Navier-Stokes Solution: Re = 10 6 

-IVI 

-N-S 




COMPRESSOR EXIT GUIDE VANE 


Comparison with Navier-Stokes Solution: Re = 10 6 

-IVI 

-N-S 



x 
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COMPRESSOR EXIT GUIDE VANE 


Loss Parameter, Q, & separation point location, 
x, ep , versus Inlet Flow Angle: 

Re = 10 6 , M_oo = 0.3 




COMPRESSOR EXIT GUIDE VANE 


Streamlines in Trailing-Edge Region: Re = 10 6 
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HIGH SPEED COMPRESSOR CASCADE 

Comparison with Navier-Stokes Solution: Re = 10 6 

-IVI 

-N-S 



HIGH SPEED COMPRESSOR CASCADE 

Comparison with Navier-Stokes Solution: Re = 10 6 

-IVI 

-N-S 



0.0 0.25 0.50 0.75 1.00 
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GOAL: UNSTEADY IVI ANALYSIS 
FOR AEROELASTIC APPLICATIONS 


• High Re unsteady cascade flows 

• Inviscid region 

- Nonlinear steady (SFLOW) => $ 

- Linearized unsteady (LINFLO) => s,v R)< £ 

note: V = + -Re{[V(^ + <j >*) + v R ] exp(iu;I)} 

- Surface conditions 

* Blades: (V - 7L) ■ n = / B {l} 

* Wakes: [V] • n = fw{$} 

• Viscous region 

- Unsteady viscous layer analysis UNSVIS 

- UNSVIS is a direct, time marching solution proce¬ 
dure 

• Inviscid/viscid coupling 

— Procedure must be developed for unsteady flows 

• Issues 

— Must modify UNSVIS to deal with moving blades 

— Matching of inviscid and viscid solutions for s and 
vr excitations 

— Need inverse unsteady viscous layer calculation 

— Inviscid/viscid coupling =► long computer run times; 
unless 

* <5 « 6 + Sexp(iu>t), i.e., linearization, or 

* Integral boundary layer calculation 
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INTERMEDIATE STEP: COUPLED 
SFLOW-IVI/LINFLO 


• Effects of strong steady interactions on unsteady 
pressure response 

• Assumptions 

- E(x, t) = S(x) + 6(x, t) 

- Strong steady inviscid/viscid interaction 

- Weak unsteady interaction 

• SFLOW-IVI will provide steady background flow infor¬ 
mation for LINFLO calculation 

- Unsteady surface pressure determined by linearized 
inviscid calculation 

- Unsteady viscous layer determined by direct 
solution procedure 
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suction 

surface 


pressure 

surface 


Wc = —0.1807, inviscid 
—0.2168, viscous 


Figure 1: LINFLO results for EGV cascade undergoing torsional vibration 
(a = (1,0), a = 0 deg, u> = 1); = 40 deg, M_= 0.30: ( ) 

inviscid; (-) viscous, Re = 10 6 . 


suction 

surface 


°f, 

Jf- 


pressure 

surface 


Wc = —0.1254, inviscid 
—0.1384, viscous 


Figure 2: LINFLO results for EGV cascade undergoing torsional vibration 

(a = (1,0), a = 0 deg, u = 1); = 54 deg, M_ 00 = 0.30: (-) 

inviscid; (-) viscous, Re = 10 6 . 






2 


suction 



Figure 1: LINFLO results for HSC cascade undergoing torsional vibration 

[a = (1,0), a = 0 deg, U) = 1]; fl_ oc = 55 deg, = 0.70: (-) 

inviscid; (-) viscous, Re = 10 6 . 



Figure 2: LINFLO results for HSC cascade undergoing torsional vibration 

[a = (1,0), a = ISO deg, u = 1]; = 55 deg, A/_<*, = 0.70: (-—) 

inviscid; (-) viscous, Re = 10 6 . 
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CONCLUDING REMARKS 


• Linearized unsteady aerodynamic analysis: LINFLO 

— Realistic 2D flow configurations 

— Arbitrary modes and frequencies of excitation 

— Efficient prediction of unsteady pressure response 

• Steady inviscid/viscid interaction analysis: SFLOW-IVI 

— 2D cascade flows 

— Local strong inviscid/viscid interactions 

— Efficient: CPU < 5 min 

— Robust: wide range of operating conditions 

• Future work 

— Transonic/supersonic gust response analysis 

— SFLOW-IVI/LINFLO coupling 

— Unsteady IVI analysis 
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N95-19382 


STEADY POTENTIAL SOLVER FOR UNSTEADY AERODYNAMIC ANALYSES 

D. Hoyniak ^ rj 

NASA Lewis Research Center 3^ ® ' 

Cleveland, Ohio 44135 

^7^3 /2- 

13 f 

Presentation Outline 

• Description of flow solver, SFLOW 

• Subsonic Calculations (Steady & Unsteady) 

Compressor Cascade (10 th Standard Configuration) 

Turbine Cascade (4 th Standard Configuration) 

GE Low Speed Research Compressor 
GE Low Speed Research Turbine 

• Transonic Calculations (Steady) 

Compressor Cascade (10™ Standard Configuration) 


Objective 

Develop steady flow solver for use with LINFLO 

• Compatible with LINFLO 

• Composite Mesh 

• Transonic Capability 


Approach 

• Steady flow potential equation written in nonconservative form 

• Newton’s Method 

• Implicit, Least-Squares, Interpolation Method used to obtain 
finite difference expressions 

• Matrix inversion routines from LINFLO 
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Differential Equations 


Steady Flow 

A 2 V 2 (J) - (y-i) V 2 ®^ - -Vcb.Zim! 

D t Dt 2 2 

= -A 2 V 2 $ + 2 

2 

= V$-V 
Dt 

Unsteady Flow 

a 2 v 2 4> - (y- 1) v 2 $>-g£ -VcD.™i =0 

£>fc DC 2 2 

-5- = io> +V$-V 
Dt 


Newton’ Method 


[Am] {<!>} = (5(0)} 
0(3f) n+1 = <J> (x) n + <|> (x) n 


Convergence Criterion 
|<(.(5c)"| <e 
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10 th Standard Configuration, Subsonic Flow Conditions 
Steady Mach Number Distribution 
M_„= 0.7, 0 .^= 55 deg 



. SFLOW 
. CASPOF 


10 th Standard Configuration, Subsonic Flow Conditions 
Unsteady Torsion Mode Response 

a. ■ 1.0, u * 0.24 a * 30 deg 
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10 th Standard Configuration, Subsonic Flow Conditions 
Unsteady Torsion Mode Response 

a * 1.0, cj * 0.24 <7 ■ 180 deg 



Standard Configuration Number 4 
Turbine Cascade Composite Mesh 
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Unsteady Pressure, Magnitude 


Standard Configuration Number 4 
Steady Surface Mach Number Distribution 



Standard Configuration Number 4 
Unsteady Aerodynamic Response 

h = (0.0016,0.0029), « - 0.24, a - -90 Deg 




Unsteady Pressure, Magnitude 


Standard Configuration Number 4 
Unsteady Aerodynamic Response 

h = (0.0016,0.0029), u - 0.24, a - 90 Deg 



GE Low Speed Research Compressor & Turbine Configurations 


Hub 


F* 


u 


I Airflow 


IGV 


Rotor 1 
R1 


nc 


Casing 


-Plane 1.0 


•Plane 1.5 


Stator 1 

81 C? 


PH 


11.43cm —\ 


Radius Hub = 0.548m Rotor 2 

Radius Casing = 0.762m 


Plane 2.0 
x Location of 
Hot-wires 
•• Location of 
Transducers 



X Location of Hot-Wires 
• • • • Location of Pressure Transducers 


Compressor Test Rig 


Turbine Test Rig 


43 

















Normalized Surface 


GE Low Speed Research Compressor 
Steady Blade Loading 


Imniisj 

Iflin 
bin 


•, ■ DATA 

I-UNSFLO 

-LEAP2D 

-UNFLO 


sisssmsi 


Chord, % 

a) Design Point (DP) 


Chord; % 

b) High Loading (HL) 


GE Low Speed Research Compressor 
Design Point, Suction Surface 


<5 0.8 [ ! 








Unsteady Pressure. Magnitude 


GE Low Speed Research Compressor 
Design Piont, Pressure Surface 



GE Low Speed Research Turbine 
Steady Blade Loading 



a) Design Point (DP) b) Low Loading (LL) 
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GE Low Speed Research Turbine 
Design Point, Pressure Surface 



o> 

c ® 

II 

U 

2 c 

co o 
CO o 

0 ) 5 

DC CO 

TJi-r 
<D C 

/A LL 


■^8 
Hi Q 




- 180 . 0 , 














Transonic Flow Calculations 


• Artificial viscosity added using rotated difference scheme 
of Jameson 

• Dissipation coefficient based on local Mach number 

• Modified Newton’s method used to solve resulting equations 


Modified Newton’ Method 
for Transonic Flow Calculations 

[AW] {<>} = {&(<!>)} 
< J >(*) n+1 = < D (*) n + ©<!>(*)'’ 


Convergence Criterion 

|<t> (*) 1 < e 


47 



10 th Standard Configuration,Transonic Flow Conditions 
M^= 0.8, £2_*=58 deg. 



10 th Standard Configuration,Transonic Flow Conditions 
Comparison with NPHASE & CASPOF Results 
M^= 0.8, Ci_ oc =58 deg. 
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Summary 

• 10 th standard configuration predictions show good agreement 
with other flow solvers 

• 4 th standard configuration turbine predictions show good 
agreement with the magnitude of measured data, however there 
are some problems with phase near trailing edge on suction 
surface 

• GE low speed research compressor and turbine predictions 
show reasonable agreement with magnitude and phase 
measurements 

• Transonic solution progressing, needs better model for artificial 
viscosity near shock, and mesh clustering capability 
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FREPS—IMPLEMENTATION 


or» 'T' 

5f 


D.V. Murthy* 

NASA Lewis Research Center 
Cleveland, Ohio 44135 


Coupling Among Disciplines 



AEROELASTIC 
ANALYSIS 




• SEQUENTIAL COUPLING 

• STRIP THEORY 

• INTERPOLATIONS 

EXTRAPOLATION 

• COORDINATE TRANSFORMATIONS 


Sequential Coupling 


GEOMETRY, OPERATING CONDITIONS 




DYNAMIC STRESSES 


•NASA Resident Research Associate at Lewis Research Center. 
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Equations of Motion 


Mechanical/Material Damping External Forcing Function from 

Forces from Damping Model Mechanical Excitation Model 





Elastic and Inertia Forces Unsteady Aerodynamic Forces 
from Structural Model from Aerodynamic Model 


Tuned System Assumption j 

CONSTANT INTERBLADE PHASE ANGLE MODES - UNCOUPLE THE EQUATIONS OF MOTION 

u,(t)=J>(t) e iM . p, = ^- r 

m n 

LEADS TO 

Mu r +Cu r +Ku r = Fj(u r (t),Ci r (t),t) +Fi(t) + Fl«(t) 

VALID FOR MOTION IN THE INTERBLADE PHASE ANGLE P r 
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Modal Transformation | 

n 

C 

^5 

M G q + C G q +Kgq =Qi(q(t),q(t),t) +0 A (t) ^Q„(t) | 


Q A (t) = Q a e lut 

Q«(t) = Q h 
q(t) = q e“ 

QX(q(t), q(t), t) = A(w) q e M 


-« 2 I + (1 m2Q [' cof ] - A(<o)] q - 0 


Discipline Communication through Databases | 


F. E. DATABASE 

GRID POINT COORDINATES, ELEMENT CONNECTIVITIES, 

NATURAL FREQUENCIES, NATURAL MODE SHAPES, NATURAL MODAL STRESSES 

STEADY AERODYNAMIC DATABASE 

AIRFOIL DEFINITION, STEADY MESH, STEADY POTENTIAL FIELD 

UNSTEADY AERODYNAMIC DATABASE 
UNSTEADY MESH 
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Input Interface 


FINITE ELEMENT ANALYSIS USER-SELECTED 

STEADY AERODYNAMIC ANALYSIS QN/ANS 

UNSTEADY AERODYNAMIC MESH GENERATION QN/ANS 

FREPS AEROELASTIC ANALYSIS KEYWORDS 

NASTRAN-LIKE 


EXAMPLES 

STRIP 1 235,236,237,238,239,240,241,242,243,244,245,246,247,248,250, 
SFLUID 100,1520.0,4635,0, 0.0000, 5618.0,1.367 
ROTOR 78 29000.0 4.790 5.290 


Output Interface 


INPUT ECHO 
GEOMETRY ECHO 
STEADY FLOW RESULTS ECHO 
FINITE ELEMENT ANALYSIS RESULTS ECHO 
UNSTEADY FLOW RESULTS FOR EACH STRIP 

GENERALIZED FORCES 

AEROELASTIC EIGENVALUES, DAMPING RATIOS, EIGENVECTORS 
FORCED DISPLACEMENT AMPLITUDES 
FORCED STRESS AMPLITUDES 
PRINTER PLOTS 
PATRAN NEUTRAL FILE 
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Current Capabilities | 

• 

INCIDENCE, CAMBER AND THICKNESS EFFECTS 

• 

SOLID/SHELL ELEMENTS 

• 

AERODYNAMIC DAMPING 

• 

DISTORTING GUST 

• 

VORTICAL, ENTROPIC, ACOUSTIC EXCITATIONS 

• 

SUBSONIC FLOW 

# 

MISTUNING SENSITIVITY MEASURE 

• 

MSC/NASTRAN INTERFACE 

• 

PRINTER PLOTS 

• 

PATRAN POST-PROCESSING 

Features in Coming Versions j 


• TIGHTER INTEGRATION OF AERODYNAMIC/AEROELASTIC MODULES 
m TRANSONIC FLOW CAPABILITY 

• DISK FLEXIBILITY / CYCLIC SYMMETRY MODEL SUPPORT 

• OTHER FEA PROGRAM SUPPORT 

• WAKE MODELING 

• VISCOUS EFFECTS (INVISCIDA/ISCID INTERACTIONS) 

• IMPROVED MISTUNING SENSITIVITY MEASURE 

• GRAPHICAL USER INTERFACE 

• 
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MEASUREMENT OF GUST RESPONSE ON A TURBINE ANNULAR CASCADE 


A.P. Kurkov and B.L. Lucci 
NASA Lewis Research Center 
Cleveland, Ohio 44135 


& 


OUTLINE 

oTest facility 

o Instrumentation 

o Wake measurements 

o Steady-state blade-surface-pressure 
measurements 

o Unsteady blade-surface-pressure measurements 
o Concluding Remarks 


INSTRUMENTATION 

UNSTEADY PRESSURE BLADE 



57 


PfWOWXNS PAGE BLANK NOT FILMED 






Frequency, engine orders 




Frequency, engine orders 
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Mean Pressure coefficient - Mean pressure coefficient 


8 



F rac tion of chord 


M = 0.27, Negative incidence 


8 
7 
6 
5 
4 
3 
2 
1 

0 
1 

-2 

0 0.2 0.4 0.6 0.8 1 

Fraction of chord 

M=0.27 Positive incidence, Reduced frequency 5 
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Suction 


0.3 

0.2 


Surface 


Port 7 





















































Frequency, 


M = 0.27, w= 10, i > 0 


Suction Surface 


0.15 

0.10 


Port 7 


















Cp : RMS Power, Far Wake, M=0.2, cu=10, i < 0 


Pressure Surface 


Suction Surface 


0.10 


Port 1 


Port 7 


















CONCLUDING REMARKS 


o At the higher Mach number, steady-state blade 
pressure distribution is to some extent dependent 
on the reduced frequency of the gust. 

o Unsteady blade pressures are strongly dependent 
on reduced frequency and incidence. Mach 
number dependence is weaker. 

o Strong pressure variation was noticed on the 
suction side of the forward part of the blade at 
lower reduced frequencies and negative 
incidence. 

o At the higher Mach number, a high-frequency 
narrow-band excitation on the suction surface 
near the trailing edge was observed. 
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EXPERIMENTAL INVESTIGATION OF AIRFOIL-GENERATED GUST FORCING FUNCTION 


S. Fleeter 
Purdue University 
West Lafayette, Indiana 47907 




Forcing Function Modeling For Flow Induced Vibration 



INTRODUCTION 


* Blade Row-Wake Interactions - Most Common Unsteady 
Aerodynamic Excitation Causing High-Cycle Blade Fatigue 



* Forced Response Design Systems Require Valid Wake Forcing 
Function Model 

* Based on Steady Performance Wake Data 



T to q « n m Dtttano. x/» 
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FORCED RESPONSE DESIGN SYSTEMS 



* Harmonic Wake - Gust Response 



1 st Harmonic Vortical Gust 

Downstream Airfoil Response 
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PROBLEM 


* Unique Requirements of Forced Response Design Systems 
* Vortical & Potential Gust Forcing Functions 



Airfoil Wake-Gust 
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RESEARCH OBJECTIVE 


* Develop Complete Model of Unsteady Aerodynamic Forcing 
Functions Within Linear Theory Framework 



TECHNICAL APPROACH 


* Perform Series of Fundamental Experiments to Investigate 
Unsteady Aerodynamic Forcing Functions 

* Vortical and Potential Gusts 
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TURBOMACHINE FLOW FIELD 



* Periodicity in Axial-Tangential Coordinate System W2 -L k 
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LINEAR THEORY GUST MODEL 


_L OB. + pV-w = 0 (Continuity) 

c§ Dt 

p^L + Vp = 0 (Momentum) 


Vortical & Potential Splitting 

W = W v + Wp 


Vortical Gust w v 


Potential Gust w p 


V-w v = 0 

DWy _ Q 

Dt 


w p = V(j) p 


P = 
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Vortical Gust Propagation 

w v = Wyexp[ -ik- (x - Wt)] 
k *W = 0 

k -Wy = 0 w v _L k 


w+ II W 

II 

ji 

w v = DWexp(-ik -x) 


• 2 Constraints 

• c|)w = 0° or 180° 

* Pw = p2 



Both 1 st 2 nd Both 

Satisfied Violated Violated Violated 

$w* 0 8 orl 80 ° P**P2 
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Potential Gust Propagation 

w p = V<j) p 

<t> p = A exp[-ik n ri + x£j 
- 7 ( 1 - 1 ^ 


Combined Vortical and Potential Gust 

W V = X D nW expf-^ktjntl + k^)] 
n 

<t>p = X A n exp[-ik nn Tl + Xn^] 



Vortical Combined Potential 


Relating Measurements to Theory 

• 2 Unknowns - A & D 

• Experimentally Measure 3 Quantities 
- u, v, p 

• Solution: 

Weighted Least Squares Approach 
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System of Equations 


X w % 

-ik^ Wyj 

Pl-W^x+w^ik^) 0 



U 

V 

P 


[T] {c} = {b} 


Weighted Least Squares 

[T] T [W] [T] {c} = [T] T [W] (b) 


where 



0 

0 

Wp 
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ROTOR 






Perforated Plate 




Measured Gust 



yortical Component Potential Component 
10*1 - 0204606 IA*I - 0.002683 




Recombine d Gust 
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Neutrally Loaded Airfoil 
Method V 




Potential Component 


4 

mbln 


ID*I ■ 0.003894 


IA*I« 0.019303 



euurany j^oaaea 
Method P 



J 


-^sL~. 


- 



0.005984 


IA*I- 0.015976 





10*1 - 0.003243 IA*1 - 0.006974 



Compressor-Loaded Airfoil 

Decayed Potential Field 

Method V 





Measured Gust Vorticai Component PfltfflTtlal ComPOnWlt 

ID*! ■ 0.005295 IA*1» 0.000675 
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Vortical Component 
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Compressor-Loaded Airfoil 
Potential Component 



0.2 0.4 0.6 0.8 1.0 # 1.2 1.4 1.0 

Axial Distance, V 


Turbine-Loaded Airfoil 
Vortical Component 
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Turbine-Loaded Airfoil 
Potential Component 



CONCLUSIONS 


* Vortical Gust Constraints Do Not Apply to Combined 
Vortical/Potential Gusts 

* Gust Should Not Be Defined Utilizing Only Unsteady Velocity Data 

* Forcing Functions Defined by Unsteady Velocity & Pressure 

* Least Squares Velocity-Pressure Method Produces Results Which 
Conform to Linear Theory 
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FORCED RESPONSE OF MISTUNED BLADED DISKS 


C. Pierre 

The University of Michigan 
Ann Arbor, Michigan 48109 


fs-oV 

y 

4 zP 


Mistuning 

■ Manufacturing tolerances, material non-uniformities, non¬ 
identical root fixtures, and in-service degradation result in 
blade-to-blade differences that destroy cyclic symmetry 


■ Small mistuning can cause large, catastrophic changes in 
blade vibrational response 

□ amplitudes of vibration of some blades may increase 
by several hundred percent, producing "rogue” blades 
and HCF failure 

□ free and forced responses may be highly sensitive to 
mistuning 

□ tuned system predictions may be qualitatively in error 
and grossly'underestimate blade forced response and 
overestimate fatigue life 


■ A credible forced response prediction system for turboma¬ 
chinery vibration must take mistuning into account 
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RESPONSE AMPLITUDE 


An Example of Mistuning Effects 
on the Free Aeroelastic Response 


□ In tuned analysis, all 
of the blades vibrate 
with equal amplitudes 
and share the excita¬ 
tion energy. 

□ Real system, with 1% 
random blade mistun¬ 
ing, has most of the 
vibration energy con¬ 
fined to 4 blades out 
of 56. 

Mistuning causes vibration localization 
—► much larger amplitudes for some blades 
—► high stresses 
—+ blade fatigue 

If unaccounted for, mistuning could cause cracks and catas¬ 
trophic blade failures. 



BLADE NUMBER 
(Pierre and Murthy, 1991) 


Effect of mistuning on forced response 
for a common blade assembly model 




Mistuning standard deviation {%) 
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Obstacles 


■ Mistimed assembly analyses are very expensive. Parametric 
studies cannot be performed 

—* need for accurate reduced-order models 

■ Mistuning is random by nature 

□ mistuning pattern (and sometimes mistuning strength) 
is typically not available 

□ mistuning differs from rotor to rotor 

□ mistuning that results from in-service degradation can¬ 
not be modeled deterministically 

—* calls for statistical and parametric tools 


■ Studies of mistuning by Afolabi, Bendiksen, Ewins, Griffin, 
Kaza, Kielb, Pierre, Sinha, Srinivasan, Mignolet, etc., have led 
to general conclusions: 

□ helps flutter 

□ increases forced response amplitudes 

■ However — quantitatively and even qualitatively different 
findings regarding other issues 

□ blade with largest amplitude 

□ forced response amplitude increase over tuned system 


■ A new perspective of the mistuning problem (Bendiksen, 
Pierre): 

□ Mistuning belongs to the broader topic of repetitive 
structures with periodicity-breaking irregularities 

□ identification of the basic physical mechanism govern¬ 
ing mistuning effects: sensitivity of aeroelastic eigen- 
solution to mistuning is inverse proportional to the dis¬ 
tance between the eigenvalues 


<5 2 A »■ 


A oj — Aofc 
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□ closeness of eigenvalues is governed by the interblade 
coupling and number of blades 

—► weakly coupled assemblies are highly sensitive 
to mistuning (interblade coupling depends on 
frequency) 

—► assemblies with many blades are more sensitive 

—> mistuning effects increase with frequency (tip 
modes) 

□ highly sensitive mistuned assemblies feature localized 
responses 


■ Formulation of a preliminary unifying theory of mistuning 


■ Demonstration of the importance of considering mistuning 
effects at the design stage 


Objectives 

Provide the designer with tools for predicting the forced 
response amplitudes of real (Le. t mistuned) bladed disks. 
Incorporate a mistuning analysis capability into forced re¬ 
sponse prediction system (FREPS) 

■ develop low-dimensional reduced-order models 

■ evaluate the significance of mistuning effects in terms 
of system parameters. Identify key parameters govern¬ 
ing sensitivity to mistuning. 

■ predict the sensitivity of the system dynamics to blade 
mistuning . 

■ determine true response amplitudes for typical mis¬ 
tuned bladed disks 

■ obtain confidence intervals for amplitudes and stresses 
and estimates of fatigue life 
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NASA research program thrusts 

■ Aeroelastic characteristics of mistimed assemblies: mode 
localization and root locus scattering 

■ Stochastic measures of sensitivity to mistuning 

□ transfer matrix based 

□ eigenvalue perturbation based 

□ localization factors 

□ composite sensitivity measure for structurally and 
aerodynamically coupled rotors 

■ Dynamics of mistuned assemblies with several component 
modes per blade. 

Effect of close blade modes on tuned and mistuned system 
dynamics. 

■ Design for low sensitivity to mistuning: formulation of an 
optimization constraint. 

■ Forced response of mistuned assemblies: 

—► physical mechanisms governing mistuning effects 
—> efficient statistical computational methods 

■ Mistuned bladed disk formulation via component mode 
analysis and validation of simple models 
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(b) Frequency Ratio ■ 0.93 



(d) Frequency Ratio ■ 0.96 



(0 Frequency Ratio ■ 0.99 













102 













Practical Significance of 
the Localization Factor 

Smid — 25 
y = 0.2 

90% amplitude decay 
by the 11th blade 
ri For y = 0.1, amplitude decays by a factor e~ 0A » 0.9 from 
one bay to the next (on average) - 

56% of the energy is transmitted to the 3rd bay 
FI Fory = 1.0, average energy transmitted to next bay is 13.5% 
and less than 0.25% of the energy reaches the 3rd bay! 

II y is an average quantity and specific realizations of mis- 
tuned systems may exhibit different behavior. 

10 

FI y can be calculated in 
terms of a universal 

Y 

sensitivity measure d,TT11d 01 
for simple models. ooi 

0.001 

FI Use in design: 

Maximum allowable localization strength y => S => corre¬ 
sponding permitted regions in parameter space. 
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Response Amplitude 


Forced Response of Mistuned Assemblies 



1900 2000 2100 2200 2300 



0.0001 0.001 0.01 0.1 1 10 100 

R 


Maximum blade amplitude in mistuned system 
Blade amplitude in tuned system 
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Closing 


■ Because of its potentially catastrophic effects such as sin¬ 
gle blade failure, mistuning must be accounted for in the 
design and analysis of blade assemblies 

■ Simple and effective mistuning capability must be incor¬ 
porated into FREPS 

■ Underlying physical mechanisms must be understood to 
generate proper reduced-order models 


Future work: 

■ Forced response: develop physical understanding and as¬ 
sociated efficient computational techniques 

■ Mistuning experiment: corroborate occurrence of local¬ 
ization and high sensitivity in nonrotating/rotating con¬ 
ditions 

■ Beneficial mistuning patterns: practical only if mistuning 
can be controlled 
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MISTUNING PATTERNS AND FORCED RESPONSE OF BLADED DISKS 


B.C. Watson 

Georgia Institute of Technology 
Atlanta, Georgia 23681 


<7) n l f 

rtf 


RESEARCH OBJECTIVES 


• BETTER UNDERSTANDING OF MISTUNED BLADED DISK 
ASSEMBLY RESPONSE CHARACTERISTICS 


• RELATIONSHIPS BETWEEN MISTUNE PATTERN, EXCITATION 
MODE, AND RESPONSE AMPLITUDE 


• OPTIMIZATION OF MISTUNING PATTERNS SUCH THAT 
FATIGUE DAMAGE IS MINIMIZED 
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BLADE/DISK MODEL 



Blade 'O' connects to blade 'N' for cyclic periodicity 


AERODYNAMIC MODEL 
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Mistune Modes 



Mean = 5 0 

1 n-i 9 

Variance = ^ Z 5f 


MODE 1 




MODE 3 
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GOVERNING EQUATIONS 


- 1 

w 

8 

td 

O 

O 

_i 


V* 

t> * 


o' 

Bjo B n B 12 


h* 

"a 

> = -i 

R, 

0 ®2i ®22 


P 

£> 


r 2 


y is the vector of disk element modal amplitudes 

* 

h Q is the vector of blade translational modal amplitudes 

cl 

a a is the vector of blade torsional modal amplitudes 


Boo _ 


-( 7 2 -7 2 fl )l + ^Y 2 h [l+[C h ]]+2Y 2 a ,[l 



B 01 =-^Y 2 h [l + [C h ]] 
B I0 =-fiJh[l + [ C h]] 

B ,, = +mYh[l+[Ch]]+Y 2 [ 4 ]] 

b i2 =y 2 [4] 
b 2 , =-y 2 K] 

B 22 =[-fl,4Y 2 I + fl,ra,Ya[ I+ [ C a]]-Y 2 ["*a]] 
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STRUCTURE OF THE CIRCULANT MISTUNE MATRIX 


[c«]= 


1 . 

—circ 

5 0 


0 


2 


5, 


-z'5 N 

1+ 2 J 





—/8 


2 + ", 

2 y 


f 


\ 


5 2 +i5^ N 
v 2+ Iy 


2 


\ 

5, +/5 N 
1 + 


SMALL MISTUNE PARAMETER 


(C h ] = e[c,,] 

[c„]=e[c„] 


Boo = Boo + e Boo 
Boi = B 0I + e B 0 , 

etc. 


By are diagonal matrices 
Bjj are circulant matrices 


ill 



RESPONSE PERTURBATION 


yI = yIo +e yIi+ e2 yl2+- 

h a — h a0 + eh a , + e h a2 H 
a a =a a0 + ea al +e 2 a a2 +-- 


PERTURBATION EQUATIONS 
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FIRST ORDER EXCITATION MODES 


FOR EXCITATION MODE p 
and 

SINGLE MISTUNE MODE r 


p, = (p-r) mod N 
p 2 = (p+r) mod N 


POSSIBLE RESPONSE MODES 

FOR MISTUNE MODE S WHICH DIVIDES N 


(p + j N/S) mod N; j = 0, 1, 2,S-l 



MODE 11 MISTUNE OFF RESONANCE 



MODE 11 MISTUNE NEAR RESONANCE 



Constant Interblade Phase Angle Mode 
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DESIGN PROBLEM 


Given a set of (mistuned) blades, find an arrangement 
(blade - slot pairing) such that the forced response 
amplitudes are minimized. 


Note: For an assembly with N blades, 
there are (N-1)! possible arrangements 


TWO PHASE SYSTEM 


• Phase One 

Continuous optimization to find Mistune Mode 
Parameters with a constraint on mistune 
strength 


• Phase Two 

Find the arrangement that best matches the 
pattern desribed by mistune mode parameters 
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CONSTRAINTS 


FOR A GIVEN SET OF BLADES, THE MEAN AND 
THE STANDARD DEVIATION ARE FIXED 


• FIXED MEAN --> DETERMINES 5 0 


• FIXED STANDARD DEVIATION ~> DETERMINES 




Mod* number 

Optimal mistune mode parameters for y = 0.99 
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Objective function values for single mode mistune cases for YM3.99 


Resoonse Amolltuda O Blade natural frequency 



1.03 



BLADE-SLOT ASSIGNMENT 


• Find {6} such that p-b\\ 2 is minimized 

• Classic Linear Sum Assignment Problem 

N-! N-I IjM 

Minimize subject to: X a > = 1 and IX = 1 

j-0 k-0 J-0 k-o 

Where: r if blade j is at slot k 

3jk 1 0; otherwise 



Best blade arrangement for optimal mistune pattern 
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0.16 



0 10 20 30 40 50 60 

Slot Number 


Best blade arrangement for mode 1 mistune pattern 
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■ Best Much 
D Alternate Pattern 

“I-1-1--1-1 

20 30 40 50 60 

Slot Number 


Best blade arrangement for alternate mistune pattern 
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Response amplitudes of best fit to mistune patterns 


LOCALIZATION 


• Localized physical blade amplitudes can imply 
non-localized modal amplitudes 


• Localized modal amplitudes implies r.on- 
localized physical blade amplitudes 
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FREPS CODE DEMONSTRATION 


M.R. Morel 
NYMA, Inc. 
Brook Park, Ohio 44142 

and 


D.V. Murthy* 

NASA Lewis Research Center 
Cleveland, Ohio 44135 


<?yY)lT 




FREPS Analysis Options 




'NASA Resident Research Associate at Lewis Research Center. 









FREPS Input/Output File Structure 
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FREPS Primary Input File 


Strip Definition 

- Aerodynamic Properties 

- Fluid and Thermodynamic Properties 

- Database for SFLOW and LINFLO Results 


jjj^Descriptior^^ottK 


■ Type of Analysis and Corresponding Parameters 
I - Flutter 

■ Flutter Frequency, Range and Iteration Parameters 

■ - Forced Response 

H Finte Element Modal Stresses; Aerodynamic/Mechanical Excitation 


Define Aerodynamic Matrix 
H and Structural Damping 


Output Request 

fjTO - Steady and Unsteady Results 
H - Line Plots 
- Raw Data 

I^^^^^ATRA^Neutra^il^^^^^ 
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FREPS Sample Input 


1 ..80 

TITLE TURBINE BLADE - FORCED RESPONSE ANALYSIS 
$ 

$ STRIP DEFINITION 
$ 

$ nodal points defining locus 


STRIP 

1 

3 

27 

75 159 243 327 411 

STREP 

1 

495 

579 

663 

747 

831 

915 

999 

STRIP 

1 

1083 

1206 

1321 

1394 

1397 

1401 

1405 

STRIP 

1 

1409 

1413 

1417 

1421 

1424 

1427 

1428 

STRIP 

1 

1426 

1422 

1419 

1414 

1410 

1406 

1402 

STRIP 

1 

1398 

1318 

1203 

1080 

996 

912 

828 

STRIP 

1 

744 

660 

576 

492 

408 

324 

240 

STRIP 

1 

156 

72 

26 

3 





FREPS Sample Input 

$ 

$ FLUID AND AERO DESCRIPTION FOR STRIPS 
$ 


$ 

fluidID 

T(°R) p(psi) 

p(lbm/in 3 ) 

a(fps) 7 

SFLUID 

100 

518.69, 14.69, 

0 . 0000 , 

1116.4, 1.4 

$ 

SROTATE 

$ 

$ 

stripID 

1 

aerolD 

e 

0.0 

W 0 

M H 


SAERO 

$ 

$ 

100 , 

0.0, 114.0, 

0.3800, 0.225 

groupID 

> stripID aeroED fluidID 


SGROUP 

100 

1 100 100 



$ 

$ STRIP DATABASE DEFINITIONS 
$ 

$ stripID filename 

DATABASE 1 airfoill 
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FREPS Sample Input 


$ 

$ ROTOR DESCRIPTION 
$ 

$ NofBLADES Sl(rpm) hubRADIUS tipRADIUS 

ROTOR 50 6000.0 4.79 5.29 

$ 

$ AEROELASTIC SET DEFINITIONS 
$ groupIDs 

AESET 100 
$ 

$ INTERBLADE PHASE ANGLE 
$ 

$ n, <r 

AESIGMA 1 72.0 

$ 

$ GUST LOADS 
$ 

$ stripED V G 

GUST 1 25.0 

$ 

$ FREQUENCY RANGE 

5 flow fhi t H Af 

FREQUENCY 9700.0 10100.0 20.00 


FREPS Sample Input 

$ 

$ NODAL DISPLACEMENTS 
$ 

$ nodes 

NODEOUT 7851 7860 

$ 

$ ELEMENT STRESSES 
$ 

$ elements 

ELEMOUT 2886 2917 2892 

$ 

$ MSC/NASTRAN PUNCH FILE OF STRESSES 
$ 

$ iounit filename 

STRFILE 1 hpl04.data 
$ 

$ FORCED RESPONSE ANALYSES REQUESTED 
$ 

RESPONSE 

u 
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j^^REP^Outpu^iles 


FREPS Output 

- Summarize Input 

- Strip Geometry & Properties 

- Steady Aerodynamic Results 

- Unsteady Aerodynamic Results 

- Aero Damping 

- Dynamic Stresses 

- Line Plots 


FREPS Raw Data 
(Optional) 

- Airfoil Contour 

- Steady Data 

- Unsteady Data 

- Root Locus Data 


^■PATRAN Neutral File 
(Optional) 

'%'M Dynamic Principal Stresses 
|fj| for the Elements 


Space Shuttle Main Engine (SSME) Blade 
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Root Locus of the Second Blade Mode (Edgewise) 
0% Modal Damping for the SSME HPOTP Blade 



1.002 


<u 1. 

j3 

re 

> 

c 

© 

o> 

iS 

«►- 

o 

to. 

re 

Q. 

re 0 , 

'5 

re 

^ 0 . 


001 


999 

998 

997 


0.996 


Root Locus of the Second Blade Mode (Edgewise) 
1% Modal Damping for the SSME HPOTP Blade 


018.46* 


25B.46* 


27.69* 

O 


36.92*. 



1—9.23* 

/ r- 304.61* 

§ O O O 


,o° 0 ^-<r 


83.08*. 

110.77* \- 350.77* 

341.54*0 / '—164.62* 

^-147.69* 


Numbers Denote 
Interblade Phase Angles 


Stable 


323.07* 


332.31* 

O 


-0.012 -0.011 -0.01 -0.009 -0.008 -0.007 

Real Part of Eigenvalue 


-0.006 
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Max Principal Stress (ksi) Max Principal Stress (ksi) 


Forced Dynamic Stresses Due to an Assumed Gust Load 
SSME HPOTP Blade 



Frequency (Hz) 


Forced Dynamic Stresses Due to Cooling Jet Excitation 
SSME HPOTP Blade 
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